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Abstract 



We give a quantum algorithm for solving instances of the satisfiability problem, based on adiabatic 
evolution. The evolution of the quantum state is governed by a time-dependent Hamiltonian that 
interpolates between an initial Hamiltonian, whose ground state is easy to construct, and a final 
Hamiltonian, whose ground state encodes the satisfying assignment. To ensure that the system 
evolves to the desired final ground state, the evolution time must be big enough. The time required 
depends on the minimum energy difference between the two lowest states of the interpolating 
Hamiltonian. We are unable to estimate this gap in general. We give some special symmetric 
cases of the satisfiability problem where the symmetry allows us to estimate the gap and we show 
that, in these cases, our algorithm runs in polynomial time. 

1 Introduction 

We present a quantum algorithm for the satisfiability problem (and other combinatorial search prob- 
lems) that works on the principle of quantum adiabatic evolution. 
An n-bit instance of satisfiability is a formula 



where each clause Ca is True or False depending on the values of some subset of the bits. For a single 
clause, involving only a few bits, it is easy to imagine constructing a quantum device that evolves 
to a state that encodes the satisfying assignments of the clause. The real difficulty, of course, lies in 
constructing a device that produces an assignment that satisfies all M clauses. 

Our algorithm is specified by an initial state in an n-qubit Hilbert space and a time-dependent 
Hamiltonian H{t) that governs the state's evolution according to the Schrodinger equation. The 
Hamiltonian takes the form 
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where each Hc\ depends only on clause Ca and acts only on the bits in Ca- H(t) is defined for t between 
and T and is slowly varying. The initial state, which is always the same and easy to construct, is 
the ground state of H{0). For each a, the ground state of Hc^{T) encodes the satisfying assignments 
of clause Ca- The ground state of H{T) encodes the satisfying assignments of the intersection of all 
the clauses. According to the adiabatic theorem, if the evolution time T is big enough, the state of the 
system at time T will be very close to the ground state of H{T), thus producing the desired solution. 
For this algorithm to be considered successful we require that T grow only polynomially in n, the 
number of bits. In this paper we analyze three examples where T grows only polynomially in n. We 
are unable to estimate the required running time T in general. 

The quantum adiabatic evolution that we are using should not be confused with cooling. For 
example, simulated annealing is a classical algorithm that attempts to find the lowest energy configu- 
ration of what we have called H{T) by generating the stochastic distribution proportional to e~^^^^\ 
where f3 is the inverse temperature, and gradually lowering the temperature to zero. In contrast, quan- 
tum adiabatic evolution forces the state of the system to remain in the ground state of the slowly 
varying H(t). 

In Section we present the building blocks of our algorithm in detail. This includes some discus- 
sion of the adiabatic theorem and level crossings. In Section || we illustrate the method on a small 
example that has three clauses, each acting on 2 bits. Each 2-bit clause has more than one satisfying 
assignment but adiabatic evolution using H{t) of the form (^^) produces the unique common sat- 
isfying assignment. In Section |^ we look at examples that grow with the number of bits in order to 
study the dependence of the required running time on the number of bits. We give three examples 
of 2-SAT problems, each of which has a regular structure, which allows us to analyze the quantum 
evolution. In these three cases the required evolution time T is only polynomially big in the number of 
bits. We also look at a version of the Grover problem that can be viewed as a relativized satisfiability 
problem. In this case our algorithm requires exponential time to produce a solution. This had to be 
so, as explained in Section 4.2. 

In Section |^ we show that our algorithm can be recast within the conventional paradigm of quantum 
computing, involving sequences of few-bit unitary operators. 



2 Adiabatic Evolution for Solving Satisfiability 

In this section we present a quantum algorithm for solving satisfiability problems. 
2.1 The Adiabatic Theorem 

A quantum system evolves according to the Schrodinger equation 

^±\i,it))^Hit)\m) (2.1) 

and the adiabatic theorem |Q| tells us how to follow this evolution in the case that H{t) is slowly 
varying. Consider a smooth one-parameter family of Hamiltonians H{s), < s < 1, and take 

H{t) = H{t/T) (2.2) 

so that T controls the rate at which H{t) varies. Define the instantaneous eigenstates and eigenvalues 
of H{s) by 



H{s) \e;s)=Ee{s) \£;s) 



(2.3) 
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with 

Eo{s) < Ei{s) < ■ ■ ■ < En-i{s) (2.4) 

where TV is the dimension of the Hilbert space. Suppose \ip{0)) is the ground state of H{0), that is, 

1^(0)) = |^ = 0;,s = 0> . (2.5) 

According to the adiabatic theorem, if the gap between the two lowest levels, Ei{s) — Eo{s), is strictly 
greater than zero for all < s < 1, then 

lim \{i = 0;s = 1\^{T))\ = 1 . (2.6) 
This means that the existence of a nonzero gap guarantees that |V'(i)) obeying (p.l[) remains very 



close to the instantaneous ground state of H{t) of the form (2.2) for all t from to T if T is big 
enough. Let us define the minimum gap by 

g^,,^ ^ mm^{Eiis) ~ Eois)) . (2.7) 
A closer look at the adiabatic theorem tells us that taking 

T » (2.8) 



9min 

where 



£ = max 

0<s<l 





dH 






ds 





(2.9) 



can make 

|(^-0;5 = 1| V(T))| (2.10) 

arbitrarily close to 1. For all of the problems that we study £ is of order a typical eigenvalue of H 
and is not too big, so the size of T is governed by 



^ mm 



2.2 The Satisfiability Problem 

Many computationally interesting problems can be recast into an equivalent problem of finding a 
variable assignment that minimizes an "energy" function. As a specific example, consider 3-SAT. An 
n-bit instance of 3-SAT is a Boolean formula, (^]^), that is specified by a collection of Boolean clauses, 
each of which involves (at most) 3 of the n bits. Each bit Zi can take the value or 1 and the i label 
runs from 1 to n. Clause C is associated with the 3 bits labeled ictJc, and kc- For each clause C we 
define an energy function 

w , r , if (z 

ic ' ^3c ' ^kc ) satisfies clause C /'o 1 1 a 

c ici jc^ kc 1^ 1 ^ if {^Zi^ , Zj^ , Zkc ) violates clause C. 

We then define the total energy h as the sum of the individual /ic's, 

h = ^hc ■ (2.12) 

c 

Clearly /i > and h{zi^ Z2, . . . , z^) = if and only if (zi, Z2, . . . , z^) satisfies all of the clauses. Thus 
finding the minimum energy configuration of h tells us if the formula has a satisfying assignment. 

We will not distinguish between conventional clauses, which compute the OR function of each 
constituent variable or negated variable, and generalized clauses, which are permitted to compute 
an arbitrary Boolean function of the constituent variables. In some of our examples it will be more 
convenient to consider generalized clauses. 
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2.3 The Problem Hamiltonian Hp 

If we go from classical to quantum computation we replace the bit Zi by a spin-^ qubit labeled by 
\zi} where Zi = 0, 1. The states \zi) arc eigenstates of the z component of the i-ih spin, 

|0)=(J) and |1)=(J) (2.13) 

so 

i(l-a«)|.,) = .,|^,) where a« = ( J _° J . (2.14) 

The Hilbert space is spanned by the N = 2'^ basis vectors \zi) \z2) ■ ■ ■ \zn)- Clause C is now associated 
with the operator Hp^c, 

Hp,c{\zi)\Z2)---\Zn)) = hc{Zic,Zjc,Zkc)\zi)\Z2)---\Zn) ■ (2.15) 

The Hamiltonian associated with all of the clauses, which we call Hp, 

Hp=Y^Hp^c (2.16) 
c 

is the sum of Hamiltonians each of which acts on a fixed number of bits. By construction. Hp is 
nonnegative, that is, Hp > for all \tjj) and Hp = if and only if is a superposition 
of states of the form l^i) \z2) ■ ■ ■ \zn) where zi, Z2, ■ ■ ■ ,Zn satisfy all of the clauses. In this context, 
solving a 3-SAT problem is equivalent to finding the ground state of a Hamiltonian. Clearly many 
other computationally interesting problems can be recast in this form. 



2.4 The Initial Hamiltonian Hb 

For a given problem, specifying Hp is straightforward but finding its ground state may be difficult. 
We now consider an n-bit Hamiltonian H^ that is also straightforward to construct but whose ground 
state is simple to find. Let H^^ be the 1-bit Hamiltonian acting on the i-th bit 

H^ = ^{l-a^^) with a« = (J J) (2.17) 

so 

H^^ \xi =x) =x\xi= x) 



where 



= 0) = J ) and |xi = 1) = i=( \ ) . (2.18) 

Continuing to take 3-SAT as our working example, clause C is associated with the bits ic, jci and 
kc- Now define 

i?B,c = iJ^*"^^ + H'i'^^ + Hi''^ (2.19) 

and 

ffB=^ffB,c. (2.20) 
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The ground state of Hb is \xi = 0) \x2 — 0) ■ ■ ■ |x„ = 0). This state, written in the z basis, is a super- 
position of ah 2" basis vectors \zi) \z2) ■ ■ ■ \z„), 

|xi=0)|x2 = 0)...|a;„ = 0) = ^^^...^|zi)|z2)---k«) • (2.21) 

Zl Z2 Z„ 

Note that we can also write 

n 

Hb^YI '^'^b (2-22) 

where di is the number of clauses in which bit i appears in the instance of 3-SAT being considered. 

The key feature of Hb is that its ground state is easy to construct. The choice we made here will 
lead to an H{t) that is of the form (|l.2| ), that is, a sum of Hamiltonians associated with each clause. 



2.5 Adiabatic Evolution 



We will now use adiabatic evolution to go from the known ground state of Hb to the unknown ground 
state of Hp. Assume for now that the ground state of Hp is unique. Consider 



so from (2.2), 



Hit) = (1 - t/T)HB + {t/T)Hp 



H{s) = (1 - s)Hb + sHp 



(2.23) 



(2.24) 



Prepare the system so that it begins at i = in the ground state of H{0) 
adiabatic theorem, if g 



Hb- According to the 



is not zero and the system evolves according to (2.1), then for T big enough 
lipiT)) will be very close to the ground state of Hp, that is, the solution to the computational problem. 

Using the explicit form of ( |2.16| ) and (|2.20| ) we see that H{t) and H{s) are sums of individual 
terms associated with each clause. For each clause C let 



Hc{t) = (1 - t/T)HBfi + {t/T)Hp,c (2.25) 

and accordingly 

Hcis) = (1 - s)Hb^c + sHp^c ■ (2.26) 

Then we have 

H{t)^Y.Hc{t) (2.27) 

c 

and 

H{s) ^Y.^c{s) . (2.28) 



This gives the explicit form of H{t) described in the Introduction as a sum of Hamiltonians associated 
with individual clauses. 
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2.6 The Size of the Minimum Gap and the 
Required Evolution Time 

Typically gmin is not zero. To see this, note from that vanishing gmm is equivalent to there being 
some value of s for which Ei{s) = Eo{s). Consider a general 2x2 Hamiltonian whose coefficients are 
functions of s 



a(s) c(s) + id{s] 

c(s) — id{s) b{s) 



(2.29) 



where a, b, c, and d are all real. The two eigenvalues of this matrix are equal for some s if and 
only if a{s) = b{s), c(s) = 0, and d{s) — 0. The curve (^a{s) , b{s) , c{s) , d{s)) in will typically not 
intersect the line 0,0) unless the Hamiltonian has special symmetry properties. For example, 



suppose the Hamiltonian (2.29) commutes with some operator, say for concreteness ax- This implies 
that a(s) = b{s) and d{s) = 0. Now for the two eigenvalues to be equal at some s we only require c to 
vanish at some s. As s varies from to 1 it would not be surprising to find c(s) cross zero so we see 
that the existence of a symmetry, that is, an operator which commutes with the Hamiltonian makes 
level crossing more commonplace. These arguments can be generalized to N x N Hamiltonians and 
we conclude that in the absence of symmetry, levels typically do not cross. We will expand on this 
point after we do some examples. 

In order for our method to be conceivably useful, it is not enough for gmin to be nonzero. We must 



be sure that gmin is not so small that the evolution time T is impractically large; see (2.8). For an 
n-bit problem we would say that adiabatic evolution can be used to solve the problem if T is less than 
for some fixed p whereas the method does not work if T is of order a" for some a > 1. Returning 



to ( |2.8| ) we see that the required running time T also depends on £ given in (2^). Using ( ^.24| ) we 



have dH/ds — Hp — Hb- Therefore £ can be no larger than the maximum eigenvalue of Hp — Hb- 



From (2.16) wc see that the spectrum of Hp is contained in {0, 1,2,..., M} where M is the number of 



terms in (2.16), that is, the number of clauses in the problem. From ( ^.22| ) we see that the spectrum 
of Hb is contained in {0, 1,2, ... ,d} where d — J^di- For 3-SAT, d is no bigger than 3M. We are 
interested in problems for which the number of clauses grows only as a polynomial in n, the number 
of bits. Thus £ grows at most like a polynomial in n and the distinction between polynomial and 
exponential running time depends entirely on ^min- 

We make no claims about the size of gmin for any problems other than the examples given in 
Section ^. We will give three examples where gmin is of order 1 /n^ so the evolution time T is polynomial 
in n. Each of these problems has a regular structure that made calculating gmin possible. However, 
the regularity of these problems also makes them classically computationally simple. The question 
of whether there are computationally difficult problems that could be solved by quantum adiabatic 
evolution we must leave to future investigation. 

2.7 The Quantum Algorithm 

We have presented a general quantum algorithm for solving SAT problems. It consists of: 



1. An easily constructible initial state (2.21), which is the ground state of Hb in (2.20) 



2. A time-dependent Hamiltonian, H{t), given by ( 2.23 ) that is easily constructible from the given 
instance of the problem; see ( 2.1(\ ) and ( ^.20 ). 



3. An evolution time T that also appears in ( 2.23| ) 



4. Schrodinger evolution according to (2.1) for time T 
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5. The final state {ipiT)) that for T big enough will be (very nearly) the ground state of Hp. 

6. A measurement of zi, Z2, • ■ • , -^n in the state \ip{T)). The result of this measurement will be 



a satisfying assignment of formula (1.1), if it has one (or more). If the formula (1.1) has no 



satisfying assignment, the result will still minimize the number of violated clauses. 

Again, the crucial question about this quantum algorithm is how big must T be in order to solve 
an interesting problem. It is not clear what the relationship is, if any, between the required size of T 
and the classical complexity of the underlying problem. The best we have been able to do is explore 
examples, which is the main subject of the rest of this paper. 



3 One-, Two-, and Three-Qubit Examples 

Here we give some one-, two-, and three-qubit examples that illustrate some of the ideas of the 
introduction. The two-qubit examples have clauses with more than one satisfying assignment and 
serve as building blocks for the three-qubit example and for the more complicated examples of the 
next section. 



3.1 One Qubit 

Consider a one-bit problem where the single clause is satisfied if and only if zi = 1. We then take 



(3.1) 



which has \zi = 1) as its ground state. For the beginning Hamiltonian we take (2.22) with ji — 1 and 
di = 1, 



(1) 



1 _ 1^(1) 

2 2 X 



(3.2) 



The smooth interpolating Hamiltonian H{s) given by ( ^.24[ ) has eigenvalues i(l ± Vl - 2s 2s'^), 
which are plotted in Fig.|l]. We see that ^min is not small and we could adiabatically evolve from 
|a;i = 0) to \zi — 1) with a modest value of T. 

At this point we can illustrate why we picked the beginning Hamiltonian, Hb, to be diagonal in 
a basis that is not the basis that diagonalizes the final problem Hamiltonian Hp . Suppose we replace 
Hb by H^ 



- 2 ~ 2^z 



(3.3) 



keeping Hp as in (^J). Now H{s) is diagonal in the z-basis for all values of s. The two eigenvalues 
are s and (1 — s), which are plotted in Fig.||. The levels cross so ^min is zero. In fact there is a 
symmetry, H{s) commutes with cr^ for all s, so the appearance of the level cross is not surprising. 
Adiabatically evolving, starting at jzi = 0), we would end up at |zi = 0), which is not the ground 
state of Hp. However, if we add to Hb any small term that is not diagonal in the z basis, we break 
the symmetry, and H{s) will have a nonzero gap for all s. For example, the Hamiltonian 



s e{l — s) 
e(l-s) l-s 



(3.4) 



has Qmin = £ for £ Small and the eigenvalues are plotted in Fig.^ for a small value of e. This "level 
repulsion" is typically seen in more complicated systems whereas level crossing is not. 
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Figure 1: The two eigenvalues of H(s) for a one-qubit example. 



0) \Z2 = 1) 

2. 



3.2 Two Qubits 

A simple two-qubit example has a single two-bit clause that allows the bit values 01 and 10 but not 00 
and 11. We call this clause "2-bit disagree." We take Hb of the form ( ^.22 ) with n — 2 and di — d2 — I, 
and we take Hp of the form ( 2.16 ) with the single 2-bit disagree clause. The instantaneous eigenvalues 
of H{s) of the form (2.24) are shown in Fig.^. There are two ground states of Hp, 
and \zi = 1) \z2 = 0). The starting state |?/'(0)), which is the ground state of H^, is (2.21) with n 
There is a bit-exchange operation |zi) \z2) \z2) \zi) that commutes with H{s). Since the starting 
state |'0(O)) is invariant under the bit-exchange operation, the state corresponding to the s = 1 end 
of the lowest level in Fig.^ is the symmetric state -^(l-Zi = 0) \z2 = 1) + \zi = 1) \z2 = 0)). The next 
level, Ei{s), begins at the antisymmetric state "^(l^^i — 0) \x2 = 1) — \xi = 1) \x2 ~ 0)) and ends at 

the antisymmetric state -^(l^i — 0) \z2 — 1) ~ \zi — 1) \z2 — 0)). Because H{s) commutes with the 
bit-exchange operation there can be no transitions from the symmetric to the antisymmetric states. 
Therefore the Ei{s) curve in Fig. ^ is irrelevant to the adiabatic evolution of the ground state and the 
relevant gap is E2{s) — Eq{s) . 

Closely related to 2-bit disagree is the "2-bit agree clause," which has 00 and 11 as satisfying 
assignments. We can obtain Hp for this problem by taking Hp for 2-bit disagree and acting with 
the operator that takes \zi) \z2) — > |zi) \z2)- Note that Hb = H^'' + H^^ is invariant under this 



transformation as is the starting state |V'(0)) given in (^.21 ). This implies that the levels of H{s) 
corresponding to 2-bit agree are the same as those for 2-bit disagree and that beginning with the 
ground state of Hb, adiabatic evolution brings you to -^(l^i — 0) \z2 = 0) + \zi = 1) \z2 = 1)). 

Another two-bit example that we will use later is the clause "imply" . Here the satisfying assign- 
ments are 00, 01, and 11. The relevant level diagram is shown in Fig.||. 
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Figure 2: The two eigenvalues of H{s) [or a one-qubit example where Hb and Hp are diagonal in the 
same basis. The levels cross so gmin — 0. 



3.3 Three Qubits 

Next we present a three-bit example that is built up from two-bit clauses so we have an instance of 
2-SAT with three bits. We take the 2-bit imply clause acting on bits 1 and 2, the 2-bit disagree clause 
acting on bits 1 and 3, and the 2-bit agree clause acting on bits 2 and 3. Although each two-bit clause 
has more than one satisfying assignment, the full problem has the unique satisfying assignment Oil. 

The corresponding quantum Hamiltonian, H{s) 
Hamiltonians each of which acts on two bits, 



(1 — s)Hb + sHp^ we write as the sum of 



Hp 



R 

imply 



H. 



13 



disagree 



Hf 



Hl^^) + {Hl,'>+Hl,^') + iH 



(3)^ 



H 



(3)n 



(3.5) 



The eigenvalues of H{s) axe shown in Fig.g. We see that (7,„i„ is not zero. Starting in the ground state 
of Hb, and evolving according to ( |2.l|) with H{t) — H{t/T) the system will end up in the ground state 
of Hp for T ^ 1/g^jj^. This example illustrates how our algorithm evolves to the unique satisfying 
assignment of several overlapping clauses even when each separate clause has more than one satisfying 
assignment. 

The alert reader may have noticed that two of the levels in Fig.O cross. This can be understood 



in terms of a symmetry. The Hamiltonian Hp of (3.5) is invariant under the unitary transformation 
V^(|zi) \z2) l^s)) — 1^2) l^i) l^a), as is Hb- Now the three states with energy equal to 4 at s = are 
\xi = 1) \x2 = 1) |a;3 ~ 0), \xi — 0) = 1) \xj, = 1), and \xi = 1) \x2 = 0) \x^ — 1). The transforma- 
tion \z) |z) in the |a;) basis is \x) (—1)^ \x) , so the states 

\xi = 1) \X2 = 1) |X3 = 0) and \xx = 0) \x2 = 1) l^g = 1) - |a;i = 1) \x2 = 0) jxa = 1) 

are invariant under V , whereas 



Cl = 0) \X2 - 1) |X3 = 1) + \XI = 1) \X2 = 0) |.T3 = 1) 
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Figure 3: A small perturbation is added to the Hamiltonian associated with Fig.^ and we see that 
the levels no longer cross. 



goes to minus itself. We call these two different transformation properties "invariant" and "odd". 
Thus at s = there are two invariant states and one odd state with energy 4. We see from Fig.^ 
that one combination of these states ends up at energy 2 when s = 1. The energy-2 state at s = 1 is 
\zi = 0) \z2 = 1) \z3 = 0), which is invariant so the level moving across from energy 4 to energy 2 is 
invariant. This means that one of the two levels that start at energy 4 and end at energy 1 is invariant 
and the other is odd. Since the Hilbert space can be decomposed into a direct sum of the invariant 
and odd subspaces and accordingly H{t) is block diagonal, the invariant and odd states are decoupled, 
and their crossing is not an unlikely occurrence. 

Since, in this simple 3-bit example, we do see levels cross you may wonder if we should expect 
to sometimes see the two lowest levels cross in more complicated examples. We now argue that we 
do not expect this to happen and even if it does occur it will not effect the evolution of the ground 
state. First note that the transformation which is a symmetry of ( |3.5| ) is not a symmetry of the 
individual terms in the sum. Thus it is unlikely that such symmetries will typically be present in more 
complicated n-bit examples. However, it is inevitable that certain instances of problems will give 
rise to Hamiltonians that are invariant under some transformation. Imagine that the transformation 
consists of bit interchange and negation (in the z basis) as in the example just discussed. Then 
the starting state |a; = 0) given by (2.21) is invariant. Assume that Hp has a unique ground state 
\zi — wi) \zi — W2) ■ ■ ■ \zn — Wn). Sincc Hp is invariant this state must transform into itself, up to a 
phase. However, from the explicit form of the ground state we see that it transforms without a phase, 
that is, it is invariant. Thus, following the evolution of the ground state we can restrict our attention 
to invariant states. The gap that matters is the smallest energy difference between the two lowest 
invariant states. 
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s 



Figure 4: The four eigenvalues of H{s) associated with "2-bit disagree". The same levels are associated 
with "2-bit agree". 

4 Examples with an Arbitrary Number of Bits 

Here we discuss four examples of n-bit instances of satisfiability. In three of the examples the problems 
are classically computationally simple to solve. These problems also have structure that we exploit 
to calculate gmin in the corresponding quantum version. In each case gmin goes like l/n^, so these 
problems can be solved in polynomial time by adiabatic quantum evolution. The other example is the 
"Grover problem" which has a single (generalized) n-bit clause with a unique satisfying assignment. 
If we assume that we treat the clause as an oracle, which may be queried but not analyzed, it takes 2" 
classical queries to find the satisfying assignment. Our quantum version has gmin of order 2~"/^, so the 
time required for quantum adiabatic evolution scales like 2", which means that there is no quantum 
speedup. Nonetheless, it is instructive to see how it is possible to evaluate (?min for the Grover problem. 

4.1 2-SAT on a Ring: Agree and Disagree 

Consider an n-bit problem with n clauses, each of which acts only on adjacent bits, that is, clause Cj 
acts on bits j and j + 1 where j runs from 1 to n and bit n + 1 is identified with bit 1 . Furthermore 
we restrict each clause to be either "agree", which means that 00 and 11 are satisfying assignments or 
"disagree", which means that 01 and 10 are satisfying assignments. Suppose there are an even number 
of disagree clauses so that a satisfying assignment on the ring exists. Clearly given the list of clauses 
it is trivial to construct the satisfying assignment. Also, if Wi,W2, ■ ■ ■ ,Wn is a satisfying assignment, 
so is wi, W2, . . . , so there are always exactly two satisfying assignments. 
The quantum version of the problem has 

Hp = H}j\ +H^l+--- + (4.1) 

where each Cj is either agree or disagree. The ground states of Hp are \wi) \w2) ■ ■ ■ \wn) and |Wi) x 
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0.9 1 



Figure 5: The four eigenvalues of H{s) associated with the 2-bit imply clause. 



1^2) • • • \wn) all in the z basis. Define the unitary transformation 

l^l) \Z2) ■■■\Zn) -> \z[) 14) •• ■ 14) 

Under this transformation Hp becomes 



z^=z 



J , if = 1 

J , if Wj = 



r23 



P afficc ' afficc 



Tjnn+l 
agree 



and the symmetric ground state of Hp is 



1 



\W) = -^(|Z1 = 0) \Z2 = 0) • • • |Z„ = 0) + |Z1 = 1) |Z2 = 1) • • • |z„ = 1)) . 



(4.2) 



(4.3) 



(4.4) 



We take Hb to be ( 2.22| ) with n bits and each di = 2. Hb is invariant under the transformation just 
given. This implies that the spectrum of H{s) = (1 — s)Hb + sHp, with Hp given by ( [f.l| ), is identical 
to the spectrum of H{s) with Hp given by (^^). Thus when we find ^min using ( |4.3| ) we will have 
found (?niin for all of the n-bit agree-disagree problems initially described. 



We can write H{s) using (4^) for Hp as 



His) = il~s) ^(1 - + i(l - a(^V(^+i)) . 



(4.5) 



We denote the s = ground state given by (2.21) as \x = 0). Define the operator G that negates the 
value of each bit in the z basis, that is, G \zi) \z2) ■ ■ ■ \zn) — |zi) [22) • • • \'Zn)- This can be written as 



(4.6) 
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0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



Figure 6: The eight levels of H{s) for the 3-bit problem with Hp and Hb given by (jg.^j. 



Since G |a; = 0) = \x = 0) and [G, H{s)\ — 0, we can restrict our attention to states that are invariant 
under G such as (4.4). 

We now write (4.5) in the invariant sector as a sum of n/2 commuting 2x2 Hamiltonians that we 
can diagonahze. First we make a standard transformation to fermion operators. To this end we define 
for j ^l,...,n, 



where 



1 -1 
2 [l -1 

It is straightforward to verify that 



A I'll 

and a+ = ^ 



where {A, B} = AB + BA. Furthermore 



for j = 1, . . . , n and 



{b]-b,)ib]^,+b,+,)^ai^^<j(^+'^ 



(4.7) 



(4.8) 
(4.9) 
(4.10) 



for j — l,...,n— 1. We need a bit more care to make sense of (4.10) for j = n. An exphcit calculation 
shows that 



(6t -6„)(6t+fei) = -Ga(")aW 



(4.11) 
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where G is given by (4.6). Since we will restrict ourselves to the G — 1 sector, ( 4.10| ) and (4.11) are 
only consistent if 6„+i = — &i, so we take this as the definition of bn+i- 



We can now reexpress H{s) of (45) in terms of the b's: 

n 

H{s) = ^{2(1 - s)blb, + ^(1 - (fot - 6^.)(&]+i + b,+i))} 



(4.12) 



Because this is invariant under the translation, bj — > and is quadratic in the bj and a 

transformation to fermion operators associated with waves running round the ring will achieve the 
desired reduction of H{s). Let 



1 " 

Pp = ^Yl e'^'P^^^bj for p = ±1, ±3, . . . , ±{n - 1) 



3 = 1 



which is equivalent to 



1 



E 



-tTTpj / n 



(4.13) 



(4.14) 



p=±l,±3,...,±(n-l) 

and is consistent with &„+i = —bi. (We assume for simplicity that n is even.) Furthermore 
and 

which follows from (^^). Substituting ( 4.14| ) into ( |4.12| ) gives 

H{s) = M^) 

p=l,3,...,(n-l) 



(4.15) 



(4.16) 



where 



+ ,s{l - cos ^ [f3l(3p - (3.p(3lp] + » sin ^ [fSl^^l - (3^(3.^] } 



(4.17) 



The Ap^s commute for different values of p so we can diagonalize each Ap separately. 

For eachp > let |ilp) be the state annihilated by both /3p and /3_p, that is, /3p \Qp) — (i-p \VLp) — 0. 
When s = 0, |ilp) is the ground state of Ap. Now Ap{s) only connects \Vtp) to jSp) = l3''_p/3jj |rip). In 
the \ flp), |Sp) basis Ap{s) is 



Ap{s) 



s + scosirp/n isisimrp/n) 
—is{simTp/n) 4 — 3s — scos7rp/n 



For each p the two eigenvalues of Ap{s) are 

i;±(s) = 2 - s ± {(2 - 3sf + 4s(l - s)(l - cos7rp/n)} ■ 



(4.18) 



(4.19) 
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The ground state energy of (4.16) is '^E~{s). The next highest energy level is E^{s) + ^ Ep{s). 

p p=3... 

The minimum gap occurs very close to s — ^ and is 



(4.20) 



for n large. 



we see that the required evolution time T must be much greater than 



Referring back to (2. 

^ Idmin ■where for this problem £ scales like n so T ^ cn^ where c is a constant. We have shown that 
for any set of agree and disagree clauses on an n-bit ring, quantum adiabatic evolution will find the 
satisfying assignment in a time which grows as a fixed power of n. 



4.2 The Grover Problem 

Here we consider the Grover problem which we recast for the present context. We have a sin- 
gle (generalized) clause, ho, which depends on all n bits with a unique (but unknown) satisfying 
assignment w = wi,W2, ■ ■ ■ , Wn- Corresponding to he is a problem Hamiltonian 

|z) , ZJ^W 

, z = w 

1 - |z = w) (z = w\ (4.21) 



Hp \z) = 



where we use the shorthand \z) — \zi) \z2) ■ ■ ■ , \zn). We imagine that we can construct H{t) = H{t/T) 
of the form ( |2.23| ) with Hp, given by ( 2.22| ) with di = \ for all i from 1 to n. Since we are evolving 
using H{t) the problem is "oracular," that is, we use no knowledge about the structure of Hp which 



could aid us in finding w other than ( [4.2l[ ). 
We can write H{s) explicitly as 

n 

H{s) = {l-s)Y^\{l-a(^^)+s{l-\z = w) {z = w\) . 
Consider the transformation given by ( |4.2| ). Under this transformation H{s) becomes 

n 

H{s) ^ (1 - s) J2 5(1 - 4'^) + -\z = 0){z = 0|) . 



(4.22) 



(4.23) 



Because the two Hamiltonians ( 4.22 ) and ( 4.25 ) are unitarily equivalent they have the same spectra 
and accordingly the same gmin- Thus it suffices to study ( 4.23| ). 

The ground state of H{0) is \x = 0), which is symmetric under the interchange of any two bits. 
Also the operator ( 4.23 ) is symmetric under the interchange of any two bits. Instead of working in the 
2"-dimensional space we can work in the (n + l)-dimensional subspace of symmetrized states. It is 
convenient (and perhaps more familiar to physicists) to define these states in terms of the total spin. 
Define S = {S^, Sy, S^) by 



Sa 



for a = x,y, z. The symmetrical states have S'^ equal to f (f + 1), where 
characterize these states as either eigenstates of Sx or Sz 



Sy - 



Sz \mz 



m 



(4.24) 
5*2. We can 

(4.25) 
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where we have suppressed the total spin label since it never changes. In terms of the z basis states 
previously introduced, 



= ^ - k) = 
for k — 0,1, . . . ,n. In particular 



l^i) 1^2) • • • \Zn) 

Z1+Z2-I hz„ = k 



^, = f ) = |z = 0) 



Now we can write H{s) in ( |4.23| ) as 

H{s) = (1 - s) [f - S.,] + s[l - \m, = f ) (m, = f |] . 



(4.26) 



(4.27) 



(4.28) 



We have reduced the problem, since H(s) is now an {n + l)-dimensional matrix whose elements we 
can simply evaluate. 
We wish to solve 



H{s)\i;)=E\i^) 



(4.29) 



for the lowest two eigenvalues at the value of s at which they are closest. Hitting ( 4.29 ) with 
(rnx = ^ — r| we get 



[s + (1 - s)r] {m^ I - r I V) ~ s {nix = f - r | = f ) (m^ = f | V') 



E{7 



We replace E by the variable A where E = s + [1 — s)\ and obtain 



Multiply by (m^ = 1 1 



and sum over r to get 



{l-s) 



n 

r — A 

r=0 



where 



Pr 



mx = r 

2 



Using (4.26) with fc = and also the identical formula with z replaced by x we have 

1 / n 



2" V r y 



(4.30) 



(4.31) 



(4.32) 



(4.33) 



(4.34) 



The eigenvalue equation (4.32) has n + 1 roots. By graphing the right-hand side of ( 4.32| ) and 
keeping < s < 1 we see that there is one root for A < 0, one root between and 1, one root between 
1 and 2, . . . , and one root between n ~ 1 and n. The two lowest eigenvalues ofi? = s + (l — s)A 
correspond to the root with A < and the root with < A < 1. We will now show that there is a 
value of s for which these two roots are both very close to zero. 
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The left-hand side of ranges over all positive values as s varies from to 1. Pick s — s* such 



that 



(1-.*) 



E 



Pr 



At s = s* the eigenvalue equation ( f4.32| ) becomes 



Pr 



A 



r{r — A) 



From (|4.34|) we know that Pq = 2"". Define u by A = 2-"/2u. Then (H.36|) becomes 



1 " 

- = E^^ 

7/ ^ ^ 



r— 1 ^ ^ 

Because of the 2~"/^ we can neglect the u piece in the denominator and we get 

1 " P 

i7->2 



which gives 



and we have 



Now 



and 



A w ± 



2fl 



(E J.2) 



-n/2 



r=l 



/ " P 

) e5 



~n/2 



So using ( |4.35|) and ( [4.40[) we have 



y^ = l + o 

r— 1 



E^ = 4 + o(4) 



2 • 2~2 



(4.35) 



(4.36) 



(4.37) 



(4.38) 



(4.39) 



(4.40) 



(4.41) 



(4.42) 



(4.43) 



which is exponentially small. 

In Fig.0 we show the two lowest eigenvalues of H{s) for the case of 12 bits. If you evolve too 
quickly the system jumps across the gap and you do not end up in the ground state of H{1). 

That grain goes like 2~"/^ means that the required time for finding the satisfying assignment grows 
like 2" and quantum adiabatic evolution is doing no better than the classical algorithm which checks 
all 2" variable assignments. In reference a Hamiltonian version of the Grover problem was studied 
with a time dependent Hamiltonian of the form 



H{t) = Hoit) + {l-\z^w){z^ w\) 



(4.44) 
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Figure 7: The two lowest eigenvalues of H{s) for the Grover problem with 12 bits. 

The goal was to choose Hjj (t) without knowing w so that Schrodinger evolution from a w-independent 
initial state would bring the system to \z = w) in time T. There it was shown how to choose Hjj so 
that the required running time T grows as 2"2 , which is then interpreted as the square-root speedup 
found by Grover. It was also shown that for any Hoit), T must be at least of order 2 "2 for the quantum 
evolution to succeed for all w. (The continuous time bound found in is closely related to the query 
bound found first in [Q.) A slight modification of the argument which gives this lower bound can be 
made for quantum evolution with 

Hit) = Hoit) + |(1 ~\z = w){z = w\) (4.45) 

and again T must be at least of order 2^ . The adiabatic evolution we studied in this section corresponds 
to Huit) = (1 — t/T)HB with Hb as described above. The lower bound just discussed shows that no 
choice of Hb can achieve better than square-root speedup. 



4.3 The Bush of Imphcations 

Ultimately we would like to know if there are general (and identifiable) features of problems which 
can tell us about the size of ^min- For the 2-SAT example of Section 4.1, gmin is of order 1/n whereas 
for the Grover problem it is of order 2~"/^. In the Grover case Hp has the property that 2" — 1 states 
have energy 1, that is, there are an exponential number of states just above the ground state. For 
the ring problem this is not so. With Hp of the form ( |4.3D there are no states with energy 1 and 
(roughly) n'^ states with energy 2. Here we present an example with an exponential number of states 
with energy 1 but for which the gap is of order l/n^. This tells us that we cannot judge the size of 
the minimum gap just from knowledge of the degeneracy of the first level above the ground state of 
Hp. 

The example we consider has n+1 bits labeled 0, 1, 2 . . . , n. There are n 2-bit imply clauses, each 
of which involves bit and one of the other n bits. Recall that the imply clause is satisfied by the bit 
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values 00, 01 and 11 but not by 10. Furthermore we have a one-bit clause that is satisfied only if bit 
has the value 1. The unique satisfying assignment of all clauses is zq = 1, zi = 1, Z2 — 1, ■ • ■ , — 1. 








1 2 3 • • • n 



Figure 8: The bush of impUcations. There is a one-bit clause that is satisGed if bit has a value of 1. 
There are n imply clauses. The j**^ imply clause is satisfied unless bit has the value 1 and bit j has 
the value 0. 

Suppose that zq = 0. Any of the 2" values of zi, Z2, ■ ■ ■ , Zn satisfy all of the imply clauses. Only 
the one bit clause is not satisfied, so these 2" variable assignments violate only one clause. There 
are n other variable assignments that violate only one clause. These have all bits set to 1 except for 
the kth bit where 1 < fc < n. In total there are 2" + n assignments that violate only one clause and 
accordingly there are an exponential number of states with energy 1. 

We can write Hp explicitly as 



,^p=l(l + ,(0))+l^(l_,(0))(l + . 



(4.46) 



To evaluate Hb from ( ^.22 ) note that bit is involved in n + 1 clauses whereas bits 1 through n are 
each involved in only one clause, so 



i?B = (n + l)i(l-a(")) + ^i(l-a«) 



i=l 



Then H{s) in terms of the spin operators (4.24) is 



i(l + af) + i(l-4°))^ + ^^ 



(4.47) 



(4.48) 



We need only consider states that are symmetrized in the bits 1 to n. We can label the relevant states 
as |zo) |mz) where zq gives the value of bit and mz labels the z component of the total spin as in 
( 4.25 ). We need to know the matrix elements of Sx in the \mz) basis. These are 



(to^ I Sx I m^) = i 



/n /n 
V2 V2 

K?( 



D-mi 



1 



ruz — 1 
^m' .m, — 1 



(4.49) 



Given (4.49) we have numerically evaluated the eigenvalues of the 2{n + l)-dimcnsional matrix 
with elements 



((4|(m;i)i/(.)(|zo)|m,)) 



(4.50) 



for values of n in the range from 20 to 120. The two lowest eigenvalues are shown in Fig.|| for n = 50. 
The gap is clearly visible. In Fig.nfl we plot log((7inin) versus log(n) and a power law dependence is 
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Figure 9: The two lowest eigenvalues of H{s) for the bush of implications with n = 50. The visible 
gap indicates that gmin is not exponentially small. 



clearly visible. We conclude that g„ 



with p 



For this problem the maximum eigenvalue 



of Hb is 2?! + 1 and the maximum eigenvalue of Hp is rt + 1, so £, which appears in at most 

grows linearly with n. Therefore we have that with T of order n^^'^'^^'' adiabatic evolution is assured. 

We also analyzed adiabatic evolution for the bush of implications using a different prescription for 
the initial Hamiltonian. We tried 



Hi, 



(4.51) 



as opposed to (2.22). This has the effect of replacing the factor of {n+ 1) in (4.47) with a 1. The effect 
on 5min is dramatic. It now appears to be exponentially small as a function of n. This means that 
with the choice of above, quantum adiabatic evolution fails to solve the bush of implications in 
polynomial time. This sensitivity to the distinction between Hb and presumably arises because 
bit is involved in (n + 1) clauses. This suggests to us that if we restrict attention to problems where 
no bit is involved in more than, say, 3 clauses, there will be no such dramatic difference between using 
Hb or H^. 



4.4 Overconstrained 2-SAT 

In this section we present another 2-SAT problem consisting entirely of agree and disagree clauses. 
This time every pair of bits is involved in a clause. We suppose the clauses are consistent, so there 



are exactly 2 satisfying assignments, as in Section 4.1. In an n-bit instance of this problem, there are 
(2) clauses, and obviously the collection of clauses is highly redundant in determining the satisfying 
assignments. We chose this example to explore whether this redundancy could lead to an extremely 
small grain- In fact, we will give numerical evidence that gmm goes like l/n^ for this problem, whose 
symmetry simplifies the analysis. 
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Figure 10: The bush of imphcations; log((7min) versus log(n) with n ranging from 20 to 120. The 
straight hne indicates that gmin ^ n~P. 



As with the problem discussed in Section 4.1, at the quantum level we can restrict our attention 
to the case of all agree clauses, and we have 



^ / J agree 
j<k 



(4.52) 



Each bit participates in [n — 1) clauses, so when constructing using (2.22) we take di — n — 1 for 
all i. We can write H{s) explicitly for this problem 



ff(.) = (l-,s)(n-l)^i(l-4^-)) 

which in terms of the total spin operators Sx and Sz is 

H{s) = (1 - s){n - 1) 



.^i(l-a«aW) 
j<k 









.2 


+ s 


— S y S y 
.4 . 



(4.53) 



(4.54) 



As in Section |4.3| , it is enough to consider the symmetric states |toz). Using (4.4£), we can find the 
matrix elements 



{m'z\H(s)\mz) 



(4.55) 



and numerically find the eigenvalues of this (n + 1) x (n + l)-dimensional matrix. 

Actually there are two ground states of i?(l), \mz = §) and \mz = corresponding to all 

bits having the value or all bits having the value 1. The Hamiltonian H{s) is invariant under the 
operation of negating all bits (in the z basis) as is the initial state given by ( ^.2l| ). Therefore we can 
restrict our attention to invariant states. In Fig.|ll]we show the two lowest invariant states for 33 bits. 
The gap is clearly visible. (i?i(0) = 64 = 2(33 — 1) because the invariant states all have an even 
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Figure 11: The two lowest eigenvalues of H{s), restricted to the invariant subspace, for overconstrained 
2-SAT with n = 33. The visible gap indicates that gmin is not exponentially small. 



number of I's in the a:-basis.) In Fig. 12 we plot log((/min) against log(n). The straight hne shows that 
5min ~ with p ~ 0.7. For this problem the maximum eigenvalues of Hb and Hp are both of order 
rt^ so £ appearing in ( ^.8|) is no larger than n^. Adiabatic evolution with T only as big as n^^^^^-* will 
succeed in finding the satisfying assignment for this set of problems. 



5 The Conventional Quantum Computing Paradigm 

The algorithm described in this paper envisages continuous-time evolution of a quantum system, gov- 
erned by a smoothly- varying time-dependent Hamiltonian. Without further development of quantum 
computing hardware, it is not clear whether this is more or less realistic than conventional quantum 
algorithms, which are described as sequences of unitary operators each acting on a small number of 
qubits. In any case, our algorithm can be recast within the conventional quantum computing paradigm 
using the technique introduced by Lloyd 

The Schrodinger equation ( ^.l| ) can be rewritten for the unitary time evolution operator U{t,to), 



and then 



i—Uit,to) = Hit)Uit,to) 



mT))^u{T,o)\m) 



(5.1) 



(5.2) 



To bring our algorithm within the conventional quantum computing paradigm we need to approximate 
U{T, 0) by a product of few-qubit unitary operators. We do this by first discrctizing the interval [0, T] 
and then applying the Trotter formula at each discrete time. 

The unitary operator U (T, 0) can be written as a product of M factors 



U{T, 0) = U{T, T - A)U{T - A, T - 2A) • • • U{A, 0) 



(5.3) 
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log(n) 

Figure 12: Overconstrained 2-SAT; log((7inin) versus log(n) with n ranging from 33 to 203. The straight 
hne indicates that gmin ~ n^. 



where A = T/M. We use the approximation 



which is vahd in (5_^) if 

\\AH{ti) 



AH{t2 



Using (2.23) this becomes 



<— for all ti,t2 G [M,(^+1)A] 



A iJp -i/B < 1 



(5.4) 



(5.5) 



(5.6) 



We previously showed (in the paragraph after Eq. ( 2.29| )) that \\Hy> — H-qW grows no faster than the 
number of clauses, which we always take to be at most polynomial in n. Thus we conclude that the 
number of factors M — T / A must be of order T times a polynomial in n. 

Each of the M terms in ( |5.3| ) we approximate as in (5^). Now H{£A) = uHb + vHp where 
u ^ 1 — {£A/T) and v — £A/T are numerical coefficients each of which is between and 1. To use 
the Trotter formula 



-iAH{eA) 



-iAuHu/K ^-iAvHr/K\K 



(5.7) 



for each£, £ = 0, 1,...,M-1, weneed A'> Af(l + A||i7B||+A||i7p||) . Since \\Hb\\ and ||Hp|| are 
at most a small multiple of the number of clauses, we see that K need not be larger than M times a 
polynomial in n. 

Now ( |5^ is a product of 2K terms each of which is e-'^^'^"^/^ or e-*^''-f^r'/^. From ( ^.22| ) we see 
that Hb is a sum of n commuting one-bit operators. Therefore e-»^"^fB/A' written (exactly) 

as a product of n one-qubit unitary operators. The operator Hp is a sum of commuting operators. 
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one for each clause. Therefore g"*^"^?/^ can be written (exactly) as a product of unitary operators, 
one for each clause acting only on the qubits involved in the clause. 

All together U{T,0) can be well approximated as a product of unitary operators each of which 
acts on a few qubits. The number of factors in the product is proportional to times a polynomial 
in n. Thus if the required T for adiabatic evolution is polynomial in n, so is the number of few-qubit 
unitary operators in the associated conventional quantum computing version of the algorithm. 

6 Outlook 

We have presented a continuous-time quantum algorithm for solving satisfiability problems, though 
we are unable to determine, in general, the required running time. The Hamiltonian that governs the 
system's evolution is constructed directly from the clauses of the formula. Each clause corresponds 
to a single term in the operator sum that is H{t). We have given several examples of special cases 
of the satisfiability problem where our algorithm runs in polynomial time. Even though these cases 
are easily seen to be classically solvable in polynomial time, our algorithm operates in an entirely 
different way from the classical one, and these examples may provide a small bit of evidence that our 
algorithm may run quickly on other, more interesting cases. 
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